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Abstract 

The Navier-Stokes equations describing laminar flow of an incom- 
pressible fluid will be solved. Different group of general solutions for 
Navier stokes equations governing Laminar incompressible fluids will 
be derived. 
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1 Introduction 

Navier stokes equations for laminar flows, in different physical conditions, 
have been solved in almost all fluid mechanics text books(l). In some phys- 
ical conditions these equations have exact solutions, but in general these 
equations are complicated nonlinear partial differential equations that can 
not be solved easily. Therefore, the authors in most cases, in order to solve 
these equations, have either used numerical methods, or they have appealed 
to Physically acceptable simplifying considerations. In this way so many ef- 
forts have been put on the recognition of physical parameters and terms that 
could be neglected without affecting the system effectively. 
Also, Navier stokes equations have a very wide application in stellar astro- 
physics, [2], [3], [4]. 
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Since 2000, The author has tried to somehow overcome the complexities gov- 
erning the analytical solutions of nonlinear partial Differential equations and 
in many respects these activities has been successful. [5], [6], [7], [8], [9], [10]. 
Since Navier stokes equations governing the incompressible fluids has a wide 
application in technology and industry, in this article we try to find explicit 
solutions of these equations for the laminar incompressible fiuids in two di- 
mensions. We have treated these set of equations in a general way without 
restricting ourselves to any especial boundary value conditions, although from 
the results found in sections 2 and 3, one can conclude that these results are 
consistent with different, logical and acceptable, boundary conditions. How 
to fit the results with a especial set of boundary conditions is beyond the 
scope of this article. We will address this issue in our next efforts. 



2 Basic equations 

Navier stokes equations in two dimensions in cartesian coordinates for lami- 
nar incompressible fluids are; 

u{x,y} ^v{x,y) ^^^n^,y}^v^ + } 



uix + vix ^ --^Pix v) + W^^^^ + 



(1) 
^) 
(2) 

du{x,y) ^ dv{x,y) ^ ^ 
dx dy 

Here, u{x, y) and v{x, y) are the x and y components of the velocity respec- 
tively, and P{x, y) stands for the pressure. The two flrst equations are the 
components of momentum equations, in the x and y direction, and obviously 
the third equation is the continuity equation. 

To solve the above set of equations, at flrst we suppose that, u and are 
functions of a float function /(x, y), that is to say, 

u = u{f{x, y)),v = v{f{x, y)) (4) 

However, this is a restricting assumption, it has no contradiction with logical 
considerations based on the following discussions. 
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1- What we do is similar to the application of separation method in solving 
the problems. At the beginning, when we apply this method we do not 
have any logical supportive argument behind that. We hope by writing the 
function in the separated form, we will come to suitable results. Of course. 
It is possible that it may not work and in fact in many cases it does not. If 
it works we have succeeded, otherwise it would be put aside. 

2- Here,What we do is to collect all the variables into a single float variable 
/. In fact, it is Some sort of variable change or transformation. 

3- Any way u and v, some how are related to each other, in advance we do not 
know this relationship, one possibility is that, they both are the function of 
the same function, this is the case, in many physical situations, for example 
in electromagnetic radiation all components of electrical and magnetic fields 
all are functions of K.r + ut and in harmonic oscillation the oscillations in 
all directions are all functions of ut .Therefore it is reasonable to think that 
the two components of the velocity are the functions of a common function 
like f{x, t). 

Now according to the above considerations equation(3) will be written as; 

du df dvdf , . 

Since / is a float and arbitrary function, so far as it fulflUs the conditions 

expressed by the equations 1- 4, any properties can be attributed to that. 
Obviously at first we prefer to chose / in such a way that the equations 
could be solved as easily as possible. We do have so Many options and as 
the first choice we assume 

(3^ = (6) 

dx dy 

with this assumption equation(5) will be simplified to 

du dv ^ 

which immediately results in 

au{f) + Pv{f) = c (8) 

As we mentioned this is not the only possible choice and for example we may 
choose / to be 

^ = (9) 
dx ^ dy ^ ' 
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or 



dx dy 

What choice is the most suitable one, it depends on the boundary conditions 
that are to be satisfied, at this point let us concentrate on equations (6) and 
(7) and look for the consequences. Prom elementary algebra The solution to 
equation(6) is , 

f = f{ax + Py) (11) 

That is to say / can be any arbitrary function of (ax+Py). From equations(7) 
and (8) we have 

dv a du d'^v a d'^u , , 
= — ^2) 

dx (5 dx'' dx^ P dx"^ 

Substituting results of (11) into equation(2), we will come to 

.du{x,y) .du{x,y) /5 1 9 d'^u{x,y) d'^u{x,y) 

(13) 

subtracting equation(l) from equation(12) will result in 

1 dP dP 

The general solution to this equation is P = P{Px — ay), that is, P can be 

any function of {fix — ay). 

Now if we rewrite equation(l) in this way; 

u[x,y) +v[x,y) v{ + J 

(15) 

it will be seen that, while according to our first assumption the left hand side 
of the above equation is a function of f{ax + f3y) the right hand side has 
to be a function of (/3x — ay), but these functions are linearly independent. 



This can not happen unless |^ = or |^ = constant. The same argument 
is true for y dependence of P; that is, ^ = or ^ = constant. In this way 
we find that our original assumption-u = u{f) and v = v{f) is only consistent 
with those physical situations, at which 

p — constant (16) 
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or 

p — (5x — ay (17) 

So far we have derived two answer for P. and each will lead to a different 
result. If we follow our analysis by taking into account the result of equation 
(15) and substituting that into the equations (1) and (2), we will see that 
both of these equations are identical. Therefore we need only one of them to 
be solved. As a result equation(l) will be written as 

du{x,y) .du{x,y) d'^u{x,y) d^u{x,y) 



since u — u{f), the following results are valid 

du df du du df du 
dx dx df dy dy df 



(19) 



dx^ df^^dx' dx^ df (9t/2 df^^dy^ dy^ df ^ ' 
From equations(6), (8), (19) and (20) will will have 

Cdud£^ ff^rf^^2 , /^^2l , , (oi) 

(5dfdy ""^df^^^dx' ^dy' ^ df^dx^ dy^^' ^ ' 

Up to this point u has been separated from other parameters. To be able to 
solve for u we have to decide about /. In the next two subsections we will 
try two different cases. 

2.1 f = ax + Py 

As we said / can be any function of ax + f3y and each / will give results that 
are consistent with an especial boundary conditions. The simplest choice will 
be 

f = ax + Py (22) 
For this choice, equation(20) will be reduced to 

c d'^u d'^u ,„ , 

(23) 



v{a^ + /?2) (i/2 df 
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The obvious solution to this differential equation is 



« = yle^^(^^""+^'^ + S (24) 

where A and B are constants of integration. Substituting the above result 
into the equation(8) v will be found out to be; 

v{f) = he -aB- aAe^^^)^'"'"^''''^) (25) 

2.2 / = e"^+^^ 
For this / we will have; 

dx dx^ ^ dy dy'^ ^ ^ ^ 

consequently equation(20) will be 

1)^ = (27) 



V(q;2 + /?2) ' df 'dp 
which has the solution 

u^A + (28) 
or in terms of x and y it will be 

u = A + B{ax + (5y) -(«"+/3^) (29) 

Accordingly from equation(8) v will be 

c — OiA 

— a(ax + (3y)^<^^W (30) 

Of course for each choice of / we will have a different solution. But the 
solutions found in this way all have the simple dependance on ax + Py. Now 
we are looking for solution that will have more complicated dependance on 
X and y.This helps us to find solutions that are consistent with other type 
of boundary conditions. In the next section we will study this aspect of the 
problem. 



6 



3 Other solutions 

In this section the basic idea is that, having two different particular solu- 
tion, how we may find other solutions. Suppose we select the two following 
particular solutions . 



u, = ^e^(^?^^"^"^"^ + B (31) 

Of course if we substitute — /3 in place of P in the above equation the result 
will also be another linear independent solution to our equations. Therefore 
our next choice will be. 

^,2^^e^^(^^""-^^^+S (32) 

We have to remind that we are quite free to chose any particular solution we 
hke. 

Now we look for those u and v and P that are functions of and u2. That 
is to say : 

U = u{Ui,U2) ,V = v{Ui,U2), P = P{Ui,U2), (33) 

Substituting equations(31) in the continuity equation, (equation(3)) will lead 
to: 

du dui du du2 dv dui dv du2 „ 

H H H = (34) 

dui dx du2 dx dui dy du2 dy 

If we differentiate Ui and U2 with respect to x and y and if we substitute the 
result in equation (32) we will have 

,^i^^^^y)^^^ + A + e^^"^-'^\a^ -P^) = (35) 

OUi OUi OU2 OU2 

Since the exponentials are linearly independent, it is obvious that we may 
write 

du ^ dv ^ , , 

c^+fi^ = ^ (37) 
au2 0U2 

Upon integration we will come to 

au + pv^ G{u2) (38) 
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au- (5v^ H{ui) (39) 

Where G and H are arbitrary functions. Prom the above equations we will 
come to 

u^^{G{u2) + H{m)) (40) 

v^^^{G{u,)-H{u,)) (41) 

As before, depending on what functions we choose for G and H, we will come 
to different results, Suppose we choose G — U2 and H — ui therefore 

u^^{ui + U2) (42) 

v = ^{ui-U2) (43) 

Substituting the above results in equations (1) and (2) we will come to 
the following results for P{x, y) 



^ P (t 1l\ I c{ax+l3y) / n c/3 ^ \ \ c{ax-(iy) i 

^ '^^ - 1/2 [-2Ae^^^^^) + [\ + e^^W^)\ {c - B)\ Ae^^(^^ ca-'iy-^ (a^ + f3^) 

(44) 

c(ax-0y) / o c/3 y \ , , _-| 

= 1/2 {-c + B)Ace^^i^^) (-1 + e ^^(^^\p-'iy-' (a'' + p"") 

(45) 



dx 



P 

If we take B — ewe will have 



= (46) 

dy 

and 

^p^-pA^e^^(^ca-'u-' [a^ + (3^)'^ (47) 

or 

Q COi X 

p = -1/2 AV^^(^^^ + (48) 

Where D is a constant of integration. 

As a result we may say that the set of functions u, v and P defined 
by equations(29),(30),(40),(41),(46) are new generations of the solutions to 
Navier stokes equations. 
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4 Discussion 



The retrieved results show that there are an infinite number of solutions to 
the Navier stokes equations. Any choice for f{ax + (3)y in equation(20) will 
result in a new solution. How to chose /, depends on the Boundary condi- 
tions to be applied. 

Since all /'s in section 2 were all specifically functions of ax + (3y the result- 
ing solutions for u and v, will also be only dependent on ax + Py. Therefore 
they will, be consistent to some especial set of boundary conditions. In sec- 
tion 3, In order to increase the capabilities of the solutions, we tried to find 
out solutions with different type of dependency on x and y. Mathematically 
,what we did in that section is somehow similar to what we do in linear 
scheme. In linear scheme we prove that if ui and U2 are particular solutions 
of a linear differential equation, then u — ciUi + C2M2 will be the general 
solution for that differential equation. In a similar manner, in section 3, we 
showed that if Ui and U2 are particular solutions to Navier stokes cqiiations 
then, u = i^{G{u2) + H{ui)) will be a general solution to that equation. 
Since from the results of section 2 there are numerous Particular solutions, 
available, therefore numerous set of general solutions are predictable for our 
set of nonlinear partial differential equations. 

The results found this way in nonlinear regime is rather different from what 
we come to, in linear discussion. In linear discussion. We argue that if, 
Ui and U2 are particular solutions, then the only general solution will be 
u — CiUi + C2U2, while here in the recent discussion we come to the result 
that for any pair of particular solutions we will have a new set of general 
solutions, consequently one can not recognise from the beginning which set 
of General solutions are consistent with the required boundary condition. 
To argue this point, is beyond the scope of this article. Here, What we may 
suggest is that; It is possible to focus from the beginning, on those float 
functions / that are consistent with all or part of the boundary conditions 
and continue our calculations based on those choices. 

In short, The scheme we followed in this article provides us numerous ac- 
ceptable solutions, and it shows that with preliminary knowledge of calculus 
analysis we can overcome the complexities of solving nonlinear set of partial 
differential equations. 

At the end, what we suggest as the topics of studies in future is; "how we 
may consider the effects of boundary conditions from the beginning" . 
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